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In a recent articled H. L. Neal presents a pedagogical ap- 
proach to density-functional theorycl in the formulation 
of Kohn and Sham,u which is still largely ignored in un- 
dergraduate teaching despite its enormous significance in 
many branches of physics, by discussing the application 
to one-dimensional two-particle systems. In this context 
Neal derives an analytic expression for the Hohenberg- 
Kohn functional F[p], given by Eq. (30) of the original 
paper, that he suggests is exact for all systems with the 
harmonic interaction u(xi, x%) = k(x\ — X2) 2 /2. The pur- 
pose of this comment is to refute this claim for arbitrary 
external potentials and to point out that the functional 
used in Ref. [I] really constitutes an approximation in the 
same spirit as the local-density approximation. 

Exploiting the universality of the Hohenberg-Kohn 
functional, Neal calculates 



F = E 



v(x)p(x) dx 



(1) 



for an analytically solvable model of two coupled har- 
monic oscillators. He then rewrites the total energy E 
and the external potential v on the right-hand side in 
terms of the density p, using exact relations that are only 
valid at the ground state of the particular model, how- 
ever. This substitution hence replaces the explicit de- 
pendence on the external potential by a system-specific 
energy surface and fails to produce a universal functional. 
In particular, the minimum of the total energy obtained 
in this way in general differs from the true ground state. 
This subtle point is best seen if the resulting functional 
F[p] is directly inserted into the variational expression 



Sp(x) 



F[p] + / (v{x) - p) p{x) dx 







(2) 



that determines the ground-state density for arbitrary 
external potentials.cl The Lagrange multiplier p enforces 
the proper normalization. In the notation of Ref. ^, which 
also provides 5F[p)/Sp(x) [Eq. (31)], one thus obtains 



p(x) = p exp 



4lo(v(x) - v{0)) 



(3) 



and hence recovers the specific relation between v and p 
that was previously employed to construct the functional 
F[p] in the first place. By design, this expression is ex- 
act for the coupled harmonic oscillators. For all other 
external potentials, however, Eq. (||) constitutes an ap- 
proximation to the ground-state density that does not 
coincide with the minimum of the true total-energy sur- 
face. 



The Kohn-Sham formalism requires F[p] as an auxil- 
iary quantity to define 



F r [p]=F[p]-T r [p], 



(4) 



where T r [p] denotes the kinetic energy of the noninteract- 
ing reference system. In the usual terminology, F r [p] con- 
tains the Hartree and the exchange-correlation part of the 
total energy, its derivative 5F r [p]/Sp(x) contributes to 
the effective potential that appears in the single-particle 
Kohn-Sham equations. It follows from the above argu- 
ment that this term is not treated exactly in Ref. [l], be- 
cause the Hohenberg-Kohn functional as well as T r [p] are 
constructed for a particular system and hence are not free 
from an implicit dependence on the harmonic-oscillator 
potential. The Hartree and exchange-correlation poten- 
tial is thus modelled by that of two coupled harmonic 
oscillators. In this sense the approach Js analogous in 
spirit to the local-density approximationEl widely used in 
atomic and solid-state physics, which similarly replaces 
the exact exchange-correlation energy by that of a homo- 
geneous electron gas with the same local density. 

Much as the local-density approximation is successful 
for weakly inhomogeneous systems, so Neal's effective po- 
tential may be applied to one-dimensional two-particle 
systems with the harmonic interaction in general confin- 
ing potentials. Indeed, Ref. ^ presents meaningful results 
for several potential wells with different shapes. For a 
quantitative analysis we carefully reexamine the exam- 
ples discussed in the original paper. We diagonalize the 
Hamiltonians using a basis of noninteracting harmonic- 
oscillator eigenfunctions. In this way all matrix elements 
can be calculated either analytically or by a numerically 
exact Gauss-Hermite quadrature. The results for case 2: 
v(x) = a\x\ and case 3: v(x) = a exp(/3x 2 ), converged 
with respect to the number of basis functions, are more 
accurate and hence differ from those quoted in Ref. [I]. 
We have set a = 1.0, (3 = 0.1 and k = 1.0. In Table | we 
contrast the exact total energy E with the Kohn-Sham 
solution £ks- Although small, the deviation is genuine. 
For comparison we also quote the estimate E var from a 
two-parameter variational wave function given in Ref. [|. 
The two approximate schemes yield similar small errors 

TABLE I. Comparison of the exact total energy E, the 
approximate Kohn-Sham solution Eks and a variational esti- 
mate -E var for two different confining potentials. 



case 


E 


-Eks 


rp a 
-C^var 


2 


1.962 81 


1.963 54 


1.967 66 


3 


2.981 49 


2.98156 


2.98168 



a From Ref. §. 
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if the quantum well resembles a harmonic potential near 
the origin (case 3), otherwise the variational wave func- 
tion is less appropriate and the effective potential gives 
better agreement with the exact solution (case 2). 

In summary, although not exact, Neal's scheme yields 
good approximate results for confined two-particle sys- 
tems with the harmonic interaction in one dimension. 
Furthermore, the procedure closely follows the spirit of 
the widely used local-density approximation and is hence 
of additional pedagogical value. We hope that it will 
find recognition and contribute towards making density- 
functional theory more accessible to students. 
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